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Abstract. We propose the formulation of a dihadron fragmentation function in 
terms of parton matrix elements. Under the collinear factorization approximation 
and facilitated by the cut- vertex technique, the two hadron inclusive cross section 
at leading order (LO) in e + e~annihilation is shown to factorize into a short 
distance parton cross section and the long distance dihadron fragmentation 
function. We also derive the DGLAP evolution equation of this function at leading 
log. The evolution equation for the non-singlet quark fragmentation function is 
solved numerically with a simple ansatz for the initial condition and results are 
presented for cases of physical interest. 
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1. Introduction 

One of the most promising signatures for the formation of a quark gluon plasma 
(QGP) in a heavy- ion collision has been that of jet quenching [Q. This phenomenon 
leads to the suppression of high px particles emanating from such collisions. Such jet 
quenching phenomena have been among the most striking experimental discoveries 
from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. 

In the investigation of jet suppression, correlations between two high px hadrons 
in azimuthal angle are used to study the change of jet structure 0- While the back- 
to-back correlations are suppressed in central Au + Au collisions, indicating parton 
energy loss, the same-side correlations remain approximately the same as in p + p 
and d + Au collisions. Given the experimental kinematics, this is considered as an 
indication of parton hadronization outside the medium. However, since the same-side 
correlation corresponds to two-hadron distribution within a single jet, the observed 
phenomenon is highly nontrivial. To answer the question as to why a parton with a 
reduced energy would give the same two-hadron distribution, one has to take a closer 
look at the single and double hadron fragmentation functions and their modification in 
medium. In this paper, we take the first step by studying the dihadron fragmentation 
function in the process of e + e~ annihilation. 

For reactions at an energy scale much above Aqcd, one can factorize the cross 
section into a short-distance parton cross section which is computable order by order 
as a scries in a s (Q 2 ) and a long-distance phenomenological object (the fragmentation 
function) which contains the non-perturbative information of parton hadronization [3] . 
These functions are process independent. If measured at one energy scale, they may 
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be predicted for all other energy scales via the Dokshitzcr-Gribov-Lipatov-Altarclli- 
Parisi (DGLAP) evolution equations In this article, we will extend this formalism 
to the double inclusive fragmentation function D^ 1,h2 (zi, z 2 ,Q 2 ) or the dihadron 
fragmentation function, where Zi,z 2 represent the forward momentum fractions of 
the two hadrons. 

2. The double fragmentation function and its evolution 

In this article, we will focus on the following semi-inclusive process 

e+ + e" -»■ 7* -> hi + h 2 + X 

of e + e _ annihilation. We consider two-jet events where both the identified hadrons h\ 
and h 2 emanate from the same jet. At leading order in the strong coupling, this occurs 
from the conversion of the virtual photon into a back-to-back quark and anti-quark 
pair which fragment into two jets of hadrons. 

In the limit of very large Q 2 of the reaction, we may invoke the collinear 
approximation. Under this approximation, at leading twist, we can demonstrate the 
factorization of the two-hadron inclusive cross section into a hard total partonic cross 
section <t and the double inclusive fragmentation function D^ lh ' 2 (z\, z 2 ) (see Ref. |H| 
for details), 

D^(z 1 ,z 2 ) + D^(z 1 ,z 2 )] . (1) 
i 

In the above equation, the leading order double inclusive fragmentation function of a 
quark is obtained as 

D^ h2 ( Zl ,z 2 ) 



In the above equation z — z\ + z 2 , p represents the momentum of the fragmenting 
parton, ph is the sum of the momenta of the two detected hadrons i.e. pi + p 2 
and n is a light-like null vector. The sum over S indicates a sum over all possible 
final hadronic states. The above equation may be represented by the diagrams of 
the cut vertex notation as that in the left panel of Fig. ^ Note that all transverse 
momentum q± up to a scale have been integrated over into the definition of the 
fragmentation function. Hadrons with transverse momenta > fi± may not emanate 
from the fragmentation of a single parton. 

In the interest of simplicity we specialize to the case of the non-singlet 
fragmentation function, D^' s 2 {zi,z 2 ) = D 1 ^ 1 '' 12 {z\, z 2 ) — D^ 1 ' 2 (zi,z 2 ). With the 
definition of the dihadron fragmentation functions in the operator formalism, which is 
shown to factorizc from the hard parton cross section at LO, we may now evaluate its 
DGLAP evolution by computing the double inclusive cross section at next to leading 
order (NLO). This constitutes a rather involved procedure; the reader is referred to 
Ref. [5] for details. The resulting DGLAP evolution consists of two parts. There 
remains the usual evolution due to the radiation of soft gluons and both detected 
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Figure 1. The left panel represents the cut- vertex representation of the dihadron 
fragmentation function. The right panel represents a Next-to-Leading order 
correction from quark and gluon single fragmentation. 



hadrons emanating with a small transverse momentum off the same parton. Hadron 
pairs with transverse momentum larger than arc produced pcrturbatively, by the 
independent single fragmentation of a quark and a gluon which emanate from the 
splitting of a quark. The contribution of this process to the evolution of the dihadron 
fragmentation function is represented by the cut-vertex diagram in the right panel of 
Fig. ^ Incorporating this diagram, we obtain the DGLAP evolution of the dihadron 
fragmentation function as 
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where the switch 1 — > 2 is meant solely for the last term. The expressions for the 
splitting functions ( P q ^ qg and P q -> qg ), single fragmentation functions, as well as the 
convolution notations (* and *) may be obtained from Ref. [S]. The above equation 
may be solved numerically and results for the cases of z\ = 2z 2 and "iz 2 are presented 
in Fig. |2 We assume a simple ansatz for the initial condition at Q 2 = 2 GeV 2 i.e. 
D^g 2 (zi, z 2 ) = D hl (zi)D h2 (z 2 ). We present results for the evolution with log(Q 2 ) at 
intervals of 1, up to log(Q 2 ) = 4.693 i.e. Q 2 = 109 GeV 2 . 



3. Discussions and Conclusions 



This study was motivated by the observation that the same side correlations of 
two high pt hadrons in central Au + Au collisions remain approximately unchanged 
as compared with that in p + p and d + Au collisions. Neglecting the differences 
in production cross section and fragmentation functions for different parton species, 



Dihadron fragmentation function 



z 1= 2z, 



12 




log(Q )=0.693 
1.693 



6 



Q 
Q 

CM 

O 



Q 2 



z 1= 3z, 



log(Q )=0.693 

1.693 

2.693 

3.693 

4.693 




Figure 2. Results of the ratio of the non-singlet quark dihadron 
fragmentation function D q % 2 (zi, Z2,Q 2 ) to the single leading fragmentation 
function Dq 1 (zi, Q 2 ). In the left panel z\ = 2z2, in the right panel z\ = 3z2- 
Results are presented from Q 2 = 2GeV 2 to 109GeV 2 . 



the integrated yield of the same side correlation over a small range of angles 
should be the ratio of the dihadron to the single hadron fragmentation functions, 
D^ lfl2 (zi, Z2, Q 2 ) / D 1 ^ 1 (z\ , Q 2 ), with Z\ and z^ being the momentum fractions of the 
triggered hadron and associated hadrons, respectively. 

Under the assumption that medium modification due to multiple scattering 
and induced gluon radiation closely resemble that of radiative corrections due to 
evolution in vacuum |7|, we evaluated the influence of DGLAP evolution on the ratio 
D^ lh2 (zi, Z2, Q 2 )/ D^ 1 (zi, Q 2 ). Our numerical results indeed show little change of the 
ratio as Q 2 is varied over a wide range of values. The evolution is shown to be strongly 
dependent, however, on the ratio of the momentum fractions of the two hadrons. In 
the results of Ref. [2] the ratio r = z\j z<i is essentially integrated over all values > 1, 
as a result the effects of the evolution are washed out. Results of the evolution of the 
more complicated singlet fragmentation functions as well as the influence of medium 
modifications will be presented in a future effort. 
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